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Abstract 

We present calculations of the spectral functions of various hadronic current correlators at fi- 
nite temperature, making use of the Nambu-Jona-Lasinio (NJL) model and the real-time finite- 
temperature formalism. We study the scalar-isoscalar correlation function in a SU(3)-flavor model, 
as well as the pseudoscalar and vector correlation functions. We relate our analysis to our recent 
calculations of the properties of mesons for T < Tc, which made use of a generalized NJL model 
that includes a covariant model of confinement. Here, we exhibit values of the spectral functions 
for a range of values of T > Tc, where it is possible to neglect the effects of confinement. We 
find important excitations in the scalar sector corresponding to what are predominately singlet 
and octet states. The singlet state, which is at 550 MeV at T = 1.2 Tc, evolves from the /o(980) 
which has an energy of about 400 MeV before it disappears from the spectrum of bound states for 
T > 0.95 Tc. The octet state seen at T = 1.2 Tc, which has a mass of about 1100 MeV, evolves 
from a nodeless state that has an energy of about 1470 MeV at T = in our model. As noted 
in the literature, these modes may play an important role in the cooling and hadronization of 
quark-gluon droplets excited in heavy-ion collisions. In the past, those researchers who have used 
the SU(2)-flavor version of the NJL (without confinement) and the imaginary-time formalism have 
found bound states for the sigma and pion for temperatures that are larger than Tc. In those works 
it is seen that the sigma and pion become degenerate after the restoration of chiral symmetry, as 
that restoration described in the NJL model. The energy of that mode then increases with increas- 
ing temperature. Our model, however, has several advantages over previous studies. We are able 
to describe meson properties below Tc, as well as the confinement-deconfinement transition. Then, 
using the real-time formalism, we are able to describe the widths of the excitations of the quark 
gluon plasma that evolve from the pion and /o(980) with increasing temperature. 

PACS numbers: 12.39.Fe, 12.38.Aw, 14.65.Bt 
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I. INTRODUCTION 



Some years ago it was suggested that the quark-gluon plasma was a "confining medium" 
in the sense that the important excitations should be color singlets [1] . At about the same 
time, Hatsuda and Kunihiro reported a study of "soft modes" of the quark-gluon plasma [2] . 
These authors used a SU(2)-flavor version of the Nambu-Jona-Lasinio (NJL) model [3] and 
considered both finite temperature and density, making use of the Matsubara imaginary- 
time formalism [4,5]. (The standard NJL model does not describe confinement, so that 
the question arises as to the modifications of the results obtained in Ref. [2] if a model of 
confinement were to be included.) In Ref. [2] it was found that the sigma meson mass, which 
lies in the continuum of the model, begins to drop in energy with increasing temperature, 
while the pion energy increases (slowly) with increasing temperature. Eventually, the sigma 
and pion become degenerate in energy, with the resulting mode increasing in energy with 
increasing temperature. These results have been confirmed in a very detailed recent study 
making use of a field-theoretic approach in the calculation of the pseudoscalar and scalar 
correlation functions at finite temperature [6]. In that work, phenomenological forms for the 
nonperturbative features of the gluon propagator were introduced and the Schwinger- Dyson 
and Bethe-Salpeter equations were solved. In many cases, the model of Ref. [6] yields results 
that arc similar to those obtained when using the NJL model. It is of interest to compare 
the results of our present study with those obtained in Ref. [6] and to note some of the 
differences between our SU(3)-fiavor analysis and the SU(2)-fiavor analysis of Ref. [6]. We 
will return to a discussion of the results obtained in our work and those obtained in Ref. [6] 
in Section VI. 

Recently, we have also seen an extensive effort in the calculation of hadronic correlation 
functions using lattice simulation of QCD at finite temperature [7-9]. These calculations 
make use of the maximum entropy method (MEM) which is reviewed in Ref. [10]. Various 
peaks are seen in the spectral functions for T < Tc which persist for T > T^. It is only at 
relatively large values of T that the correlation functions go over to the smooth behavior 
expected for a weakly interacting system [7-9] . 

Although we have introduced a critical temperature for the purposes of our discussion, 
we should note that the quark masses and vacuum condensates do not go to precisely zero at 
large temperature, since the current quark mass is always present in the model. Indeed, the 



3 




FIG. 1: Temperature-dependent constituent mass values, mu{T) and ms{T), calculated in a mean- 
field approximation [12] are shown. [See Eq. (2.2)]. Here m° = 0.0055 GeV, m° = 0.120 GeV, and 
G{T) = 5.691[1 - 0.17(r/rc)] GeV^, if we use Klevansky's notation [12]. (The value of G used in 
our work is defined such that it is twice the value of G used in Ref. [12].) 

presence of a small current mass, = 0.0055 GeV, makes for a significant change in the 
behavior of m„(T) with increasing temperature. On the other hand, it is useful to introduce 
a value of of the order of 170 MeV to facihtate the discussion. The significance of that 

value for the characterization of the constituent mass of the up (or down) quark may be 
inferred from Fig. 1, which exhibits results calculated for m° = 0.0055 GeV and m° = 0.120 
GeV. 

It is of interest to obtain insight into the lattice results for spectral functions by study- 
ing hadronic current correlators in a model, such as the NJL model, that describes the 
restoration of chiral symmetry with increasing temperature. We present such a study in 
this work. The organization of our work is as follows. In Section II we review our model 
of confinement and our calculations of meson properties for T < Tg. In Section III we de- 
scribe the calculation of the imaginary part of the vacuum polarization function of the NJL 
model, when use is made of the real-time finite-temperature formalism. Once the imaginary 
part is calculated, the real part may be obtained using a dispersion relation. In Section 
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IV we show how hadronic current correlation functions may be constructed in terms of the 
vacuum polarization functions calculated in Section III. In Section V we present results of 
numerical calculations of the various hadronic current correlation functions considered in 
this work. Finally, Section VI contains some further discussion and some comparison of the 
Minkowski-space analysis of the present work and the Euclidean-space analysis of Ref . [6] . 

II. CALCULATION OF MESON PROPERTIES AT FINITE TEMPERATURE IN 
A GENERALIZED NJL MODEL WITH CONFINEMENT 

It is useful to record the Lagrangian used in our calculations of meson properties 

i=0 

+^{det[^(l + 75)?] + det[^(l - 75)?]} 

+^conf , (2.1) 

Here, mP is a current quark mass matrix, mP — diag(m^, m^, m°). The A' are the Gell- 
Mann (flavor) matrices, A° = a/2/3 1 with 1 being the unit matrix. The fourth term on the 
right-hand side of Eq. (2.1) is the 't Hooft interaction. Finally, Cconf represents the model 
of confinement we have used in our work. 

As noted earlier, we have recently reported results of our calculations of the temperature 
dependence of the spectra of various mesons [11]. These calculations were made using our 
generalized NJL model which includes a covariant model of confinement. We have presented 
results for the tt, K, ao, /o and mesons in Ref. [11]. 

There arc three important temperature-dependent features of our model. Temperature- 
dependent constituent quark masses were calculated using the equation [12] 

m(T) P2 1 

m(T) =m' + 2Gs{T)N,^-l + / dp— tanh(-/3Ep) . (2.2) 

Jo ^ 

Here, mP is the current quark mass, Gs{T) is a temperature-dependent coupling constant 
introduced in our model [11] 



Gs{T) = Gs 



1-0.17|| 



(2.3) 
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A^c = 3 is the number of colors, (3 = l/T and Ep = + m^(T)]^^^ . Further, A = 0.631 
GeV is a cutoff such that |p| < A. Results obtained for the up (or down) and strange quark 
masses are given in Fig. 1. In calculating the constituent mass values we have neglected the 
confining interaction. That interaction was taken into account in our earlier Euclidean- space 
calculation of the quark self-energy [13], which also included the effects related to 't Hooft 
interaction. We found that, to a good approximation, we could neglect the confining and 't 
Hooft interactions, if we modified the value of the NJL coupling constant, Gs, and we adopt 
that approach when using Eq. (2.2). 

In addition to the temperature dependence of the coupling constant and constituent 
mass values, we also introduced a temperature-dependent confining potential, whose form 
was motivated by recent lattice simulations of QCD in which the temperature dependence of 
the confining interaction was calculated with dynamical quarks [14] . (See Fig. 2.) In order to 
include such effects, we modified the form of our confining interaction, V'-^{r) — «;r exp[— //r], 
by replacing /i by 

t^iT) = T ^\ ..1 > (2-4) 



1-0.7 



with Ho — 0.010 GeV. The maximum value of V'^{r, T) is then 

M = (2.5) 

^ ^[1-0.7(T/Te)^] ^ 
/ioe 

with Tmax — ^ / f^{T) ■ To better represent the qualitative features of the results shown in Fig. 
2, we can use V^{r, T) = rer exp[-//(r)r] for r < r^ax and V^{r, T) = V^^{T) for r > r^ax- 
We also note that we use Lorentz- vector confinement and carry out all our calculations in 
momentum-space. The value of k used in our work is 0.055 GeV^. 

We do not attempt to review the details of our calculations of temperature-dependent 
meson spectra [11]. However, in Figs. 3 and 4 we show the results for the tt and /o mesons 
obtained in Ref. [11]. (These calculations were made with the temperature-dependent 
potentials shown in Fig. 5.) In Figs. 3 and 4, it may be seen that, as T ^ Tc, fewer states 
are bound, with no bound states remaining at T = Tc ■ Our confining potential is finite 
aXT — Tc, but has a zero value for T — 1.195Tc. (The potential is defined to be equal 
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FIG. 2: A comparison of quenched (open symbols) and unquenched (filled symbols) results for the 
interquark potential at finite temperature [14] . The dotted line is the zero temperature quenched 
potential. Here, the symbols for T = O.SOTc [open triangle], T = 0.88Tc [open circle], T = O.SOTc 
[open square], represent the quenched results. The results with dynamical fermions are given 
at T = 0.68Tc [solid downward-pointing triangle], T = O.SOTc [solid upward-pointing triangle], 
T = 0.88Tc [solid circle], and T = OMTc [solid square]. 

to zero for T > l.lQSTc-) However, it is important to note that Gs{T) — only for the 
quite large value T — 5.88Tc . (The couphng constants of the model are taken to be zero for 
T > 5.88Tc.) Even though the confining potential is absent for T > 1.2Tc, the continued 
presence of the short-range NJL interaction has important consequences for the dynamics 
of the quark-gluon plasma, as we will explore in the following sections. 



III. POLARIZATION FUNCTIONS AT FINITE TEMPERATURE 

The basic polarization functions that are calculated in the NJL model are shown in 

— * 

Fig. 6. We will consider calculations of such functions in the frame where P = 0. In 
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FIG. 3: The mass values of the pionic states calculated in Ref. [11] with Gt^{T) = 13.49[1 — 

0. 1TT/Tc] GeV-2, Gy(T) = 11.46[1 - O.ITT/Tc] GeV-^, and the quark mass values given in Fig. 

1. The value of the pion mass is 0.223 GeV at T/T^ = 0.90, where m„(T) = 0.102 GeV and 
ms{T) = 0.449 GeV. 

our earlier work, calculations were made after a confinement vertex was included. That 
vertex is represented by the filled triangular region in Fig. 6. However, we here consider 
calculations for T > 1.2Tc where confinement may be neglected. We will, however, use the 
temperature-dependent mass values shown in Fig. 1. 

The procedure we adopt is based upon the real-time finite-temperature formalism, in 
which the imaginary part of the polarization function may be calculated. Then, the real 
part of the function is obtained using a dispersion relation. The result we need for this 
work has been already given in the work of Kobes and Semenoff [15]. (In Ref. [15] the 
quark momentum in Fig. 6 is /c and the antiquark momentum is k — P. We will adopt that 
notation in this section for ease of reference to the results presented in Ref. [15].) With 
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FIG. 4: Mass values of the /o mesons calculated in Ref. [11] with Goo{T) = 14.25[1 - O.ITT/Tc] 
GeV-2, GssiT) = 10.65[l-0.17r/rc] GeV-^, GosiT) = 0.495[l-0.17r/rc] GeV-^, and G8o(r) = 
Gos{T) in a singlet-octet representation. The quark mass values used are shown in Fig. 1. The fo 
has a mass of 0.400 GeV at T/Tc = 0.95. 

reference to Eq. (5.4) of Ref. [15], we write the imaginary part of the scalar polarization 
function as 

{(1 - n,{k) - n^ikMP^ - E^{k) - E^ik)) 
-iMk) - n2(A;))5(P° + E,{k) - E^ik)) 
-{n^ik) - ni{k))5iP° - E,{k) + E^ik)) 
-(1 - ni{k) - n2{k))d{P'' + Ei{k) + E2{k))} . 

Here, Ei{k) = [k^ + m?(r)]V2. Relative to Eq. (5.4) of Ref. [15], we have changed 
the sign, removed a factor of g'^ and have included a statistical factor of 2Nc, where the 
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FIG. 5: The potential V'^{r,T) is shown for T/Tc = [sohd hnc], T/Tc = 0.4 [dotted hne], T/Tc = 
0.6 [dashed Una], T/Tc = 0.8 [dashed-dotted hne], T/Tc = 0.9 [short dashes], T/Tc = 1-0 [dashed- 
(double) dotted hne]. Here, V'^{r,T) = exp[-/x(r)r], with /x(r) = 0.01GeV/[l - (}.7{T/Tcf] 
and K = 0.055 GeV^. 

factor of 2 arises from the flavor trace. In addition, we have included a Gaussian regulator, 
exp[— ^^/cK^], with a — 0.605 GeV, which is the same as that used in most of our apphcations 
of the NJL model in the calculation of meson properties [16-21]. We also note that 

and 

For the calculation of the imaginary part of the polarization function, we may put A;^ = 
m1{T) and {k — P)^ = m^iT), since in that calculation the quark and antiquark are on- 
mass-shell. In Eq. (3.1) the factor Ps arises from a trace involving Dirac matrices, such 
that 

(3s = -Tr[{lt + mi)(^ - f + m2)] (3.4) 
= 2p2-2(ml + m2)^ (3.5) 

10 



hJ(p') 




P- 



p... 



o 
o 



...p 



...p 



FIG. 6: The upper figure represents the basic polarization diagram of the NJL model in which 
the lines represent a constituent quark and a constituent antiquark. The lower figure shows a 
confinement vertex [filled triangular region] used in our earlier work. For the present work we 
neglect confinement for T > 1.2Tc, with = 170 MeV. 

where nii and 1112 depend upon temperature. In the frame where P = 0, and in the case 
mi = 7712, we have = 2Po (1 ~ Am^ / Pq). For the scalar case, with 77ii = 7712, we find 



We may evaluate Eq. (3.6) for m{T) = 77i„(r) = md{T) and define Im Jf (P^T). When 
we put 777(T) = ms{T), we define Im Jg{P'^, T). These two functions will be needed for our 
calculation of the scalar-isoscalar correlator discussed in the next section. The real parts 
of the functions J^(P^,T) and Jf(P^,T) may be obtained using a dispersion relation, as 
noted earlier. 

For pseudoscalar mesons, we replace (5s by 




(3.6) 



where 
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(3.7) 




(3.8) 
(3.9) 



2P' - 2(7711 - m-2f , 



which for 777i = 7772 is (5p — 2Pq in the frame where P 



0. We find, for the tt mesons. 




2\ 1/2 



Im J^(P2,T) 



e 



/"'[I - 2771 (A;)], 



(3.10) 
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where k"^ — Pq/4: — m^(T), as above. Thus, we see that, relative to the scalar case, the 
phase space factor has an exponent of 1/2 corresponding to a s-wave amplitude. For the 
scalars, the exponent of the phase-space factor is 3/2, as seen in Eq. (3.6). 
For a study of vector mesons we consider 

P^u = Trb^d^ + mi)7.(^ -r + ^2)] , (3.11) 

and calculate 

g'^'P^^^^P^ -ml-ml + 4:mim2], (3.12) 

which, in the equal- mass case, is equal to APq + 8m^(r), when mi = 1712 and P = 0. This 
result will be needed when we calculate the correlator of vector currents in the next section. 
Note that for the elevated temperatures considered in this work mu{T) — md{T) is quite 
small, so that AP^ -\- 8m1{T) can be approximated by APq when we consider the p meson. 



IV. CALCULATION OF HADRONIC CURRENT CORRELATION FUNCTIONS 

In this section we consider the calculation of temperature-dependent hadronic current 
correlation functions. The general form of the correlator is a transform of a time-ordered 
product of currents, 

^C(P^^) = j d'xe'''-^ « T{j{x)m » , (4.1) 

where the double bracket is a reminder that wc arc considering the finite temperature case. 

For the study of pscudoscalar states, we may consider currents of the form jp,j(x) = 
q{x)ij5X^q{x), where, in the case of the tt mesons, i — 1,2, and 3. For the study of scalar- 
isoscalar mesons, we introduce js,i{x) — q{x)X'q{x), where i = for the flavor-singlet current 
and i = 8 for the flavor-octet current. 

In the case of the tt mesons, the correlator may be expressed in terms of the basic vacuum 
polarization function of the NJL model, Jp{P^,T) [12, 22, 23]. Thus, 

where Gt^{T) is the coupling constant appropriate for our study of the tt mesons. We have 
found G7r(0) = 13.51 GeV~^ by fitting the pion mass in a calculation made at T = 0, with 
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nT'u — TTT'd — 0.364 GeV [11]. (We will describe the calculation of Jp{P'^,T) later in this 
section.) 

The calculation of the correlator for scalar-isoscalar states is more complex, since there are 
both fiavor-singlet and flavor-octet states to consider. We may define polarization functions 
for u, d and s quarks: Ju{P'^,T), Jd{P'^,T) and Js{P^,T). (These functions do not contain 
the factor of 2 that would otherwise appear when forming the flavor trace.) In terms of 
these polarization functions we may then define 



Joo(P^T) = -[J„(P^^) + MP',T) + Js{P\T)] , 



(4.3) 



MP',T) = ^[J^{P\T) + Jd{P\T) - 2Js{P\T)] , 



(4.4) 



and 



Jss{P', T) = -[Ju{P\ T) + UP\ T) + AUP\ T)] . 



We also introduce the matrices 

J(p^T) 



Joo(P^^) Jo8(P^T) 

J8o(p^T) J88(p^T) 



(4.5) 



(4.6) 



G(T) = 



G'oo(^) GosiT) 
G8o(T) G88(T) 



and 



C{P^,T) 



We then write the matrix relation 



Coo(P^T) Co8(P^T) 
CsoiP^T) CssiP^T) 



(4.7) 



(4.8) 



C(P^ T) = J(P^ T) [1 - G(T)J(P^, T)]-^ . 



(4.9) 



For some purposes it may be useful to also define a t matrix 



t(P',T) = [1 - G{T)J{P',T)]-'G{T) , 



(4.10) 



where t{P'^,T) has the structure shown in Eqs. (4.6)-(4.8). The same resonant structures 
are seen in both C(P^ T) and i(p2, T). 
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For a study of the correlators related to the p meson, we introduce conserved vector 
currents jfjL,i{x) — q{x)^ij_\iq{x) with i — 1,2 and 3. In this case we define 

J^^{P\T) ^[g^""- J,{P\T) (4.11) 



and 



C;''{P\T) = (^g^'' - C,{P\T) , (4.12) 

taking into account the fact that the current j^,i{x) is conserved. We may then use the fact 
that 

J,(P^T) = \g,,j;'^{P\T) (4.13) 



2Nr 



+ 2ml{T) 



3 

2 
3 



2(rr\ \ 1/2 

H-///, 

1 - 



1^) e-'y''[l-2n,(k)] (4.14) 



47r 

UP\T). (4.15) 



See Eq. (3.7) for the specification oi k — \k\. We then have 

V. RESULTS OF NUMERICAL CALCULATIONS 

We begin our presentation with the results we have obtained for scalar-isoscalar excita- 
tions. In Fig. 7 we show the values of Im Joo(P^) for T/T^ = 1.2, 1.6, 2.0, 4.0 and 6.0. 
Figure 8 exhibits similar results for Im Jo8(P^), while Fig. 9 exhibits the results obtained 
for Im J88(-P^)- The real parts of these functions arc obtained using a dispersion relation. 
Once we have the real and imaginary parts of these functions, we can calculate the elements 
of the correlation functions. 

In Figs. 10-12 we show Im Coo(P^), Im Cq^^P"^) and Im C^s^P'^). In the case of Im Coo(P^) 
we see a sharp peak at 550 MeV. This peak represents the lowest /o state at that temperature 
and has as its "parent" the /o(980) meson. [See Fig. 4] The decay of the state at 550 MeV 
to two pions will be suppressed since the pionic excitation see in Fig. 13 is essentially 
degenerate with the /q excitation for T > T^. (Note that in the fiavor-SU(3) version of 
the NJL model the /o(980) replaces the sigma of the fiavor-SU(2) NJL model as the chiral 
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FIG. 7: Values of Im Joo(P^) are shown. Here, T/T^ = 1.2 [solid line], 1.6 [dashed line], 2.0 [dotted 
line], 4.0 [dashed-dotted line] and 6.0 [dashed- (double) dotted line]. 
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FIG. 8: Values of Im Jo8(P^) are shown. [See the caption of Fig. 7.] 

partner of the pion. Thus, it is the state that evolves from the /o(980) with increasing 
temperature that becomes degenerate with the pion upon restoration of chiral symmetry.) 

The peak seen at about 1100 MeV in Fig. 12 is mainly an octet state and is quite 
prominent in ImC88(-P^)- When we calculate the mixing angle for the state at 550 MeV, 
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FIG. 10: Values of lmCoo(^'^) are shown. [See the caption of Fig. 7.] 

we find 6i = 27.7° corresponding to the state being about 78% flavor- singlet. A calculation 
of the mixing angle for the state at 1100 MeV yields 62 = 17.8° which makes the state 95% 
flavor-octet. These mixing angles appear in the schematic representation 

|/o(247)) = cos^iAo + sinMs, (5.1) 
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FIG. 11: Values of IuiCqs^P^) are shown. [See the caption of Fig. 7.] 



0.12 - 




P' (GeV') 



FIG. 12: Values of ImCssiP^) shown. [See the caption of Fig. 7.] 

|/o(860)) = - sin 62X0 + cos 62X8 • (5.2) 

It may be seen in Fig. 10 that the lowest /o state moves up in energy and becomes wider 
when T = l.GTc, while the state that is at 1100 McV at T = 1.2Tc moves down in energy 
at T = l.QTc. (See Fig. 12.) We also see that as the temperature is increased, the states 
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FIG. 13: Values of ImC7r(-P^) axe shown. [See the caption of Fig. 7.] 

broaden further and the curves eventually become rather featureless. 

In Fig. 13 we show Im C^(P^) for the same set of temperatures used to exhibit the 
properties of the scalar-isoscalar excitations. The curve for T — 1.2Tc [solid hne] in Fig. 13 
has a peak at 547 MeV. With increasing temperature the excitation becomes much broader. 

Results for Im Cp(P^) are shown in Fig. 14. There is a peak at about 700 MeV that has 
as a "parent" the p(770). Again, we see increased broadening with increasing values of the 
temperature. 

VI. DISCUSSION 

The formalism we have developed allows for a description of the mass values of various 
mesons and their radial excitations for the range of temperatures T < T^. Our model 
also describes the confinement-deconfinement transition at T = T^. When we calculated 
hadronic current correlators in the real-time formalism we are able to describe the widths 
of the excitations for T > Tc. It is of interest to contrast our results with those obtained 
using the imaginary-time formahsm [4, 6, 22]. The results for the behavior of the sigma and 
pion masses are similar in the NJL model [4] and the field-theoretic model of Ref. [6]. In 
Ref. [6] the value of Tc ~ 150 MeV is obtained when the current quark mass is zero. For a 
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FIG. 14: Values of Im Cp(P^) are shown. [See the caption of Fig. 7.] 

finite value of the current quark mass, we may make reference to Fig. 6 of Ref . [6] . There, 
the sigma mass has a minimum value of approximately 255 MeV at T = 160 MeV. At that 
temperature the pion mass is about 220 MeV. If we consider T — 1.2Tc = 180 MeV, we 
have m^r = m^. = 440 MeV in Ref. [6]. In our work, we have both mj„ and at about 
550 MeV when T = 1.2Tc. Probably, one of the more significant differences between the 
imaginary-time formalism used in Refs. [4] and [6] and the real-time formalism used here is 
the absence of information concerning the widths of the excitations in the imaginary-time 
formalism. As may be seen in our Figs. 7-14, the widths of the excitations may play a very 
important role in understanding the properties of the quark-gluon plasma. 

It is beheved that measurements of particle ratios obtained in central heavy-ion coUisions 
contain useful information concerning the properties of the quark-gluon plasma. In a recent 
study, different assumptions concerning the dynamics of the plasma were made and their 
influence upon predicted particle ratios were obtained [24] . The authors consider calculations 
based upon both a noninteracting gas model and a chiral SU(3) a — uj model. They find 
that "the extracted chemical freeze-out parameters differ considerably from those obtained 
in simple noninteracting gas calculations..." when the a — ui model is used to calculate 
particle ratios. We believe our model may be useful in performing calculations of the type 
reported in Ref. [24]. 
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